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Numerical calculations are presented, which relate the intensity and polarization of Jupiter's 
decimeter radiation to the distribution of synchrotron-radiating electrons in its "Van Allen belts." 
The calculations are based on the simple model of a dipole magnetic field centered in the planet but 
inclined to the axis of rotation. Shadowing by the planet's disk is taken into account. An appropriate 
choice of the parameters of the model enables one to account for (1) the intensity, (2) the spectrum, 
(3) the beaming, and (4) the degree of polarization of Jupiter's decimeter radiation. However, the model 
cannot account for the observed asymmetries in the beaming and polarization. 



1. Introduction 

Nearly 4 years ago the author undertook a study of 
the theory of synchrotron radiation from stars with 
dipole magnetic fields (results reported in [Thorne, 
1963]). As part of this study a model was constructed 
in which the beaming of the radiation was determined 
by the pitch-angle distribution and energy distribution 
of the radiating electrons. Only after this study was 
completed did the author become aware of the dis- 
covery by Morris and Berge [1962] of the obliquity of 
Jupiter's magnetic field and of the consequent ability 
of Earth-based observers to study the beaming of 
Jupiter's decimeter radiation. However, J. A. Roberts 
and M. M. KomesarofT realized immediately that the 
author's calculations for stars are also applicable 
to Jupiter; and they communicated with the author 
concerning the interpretation of the decimeter radia- 
tion in the light of those calculations. In order to make 
the calculations more applicable to Jupiter, the author 
modified the original stellar model to take into account 
shadowing by the planet's disk and the absence of 
electrons whose mirror points would be inside the 
planet. The results of calculations based on this 
modified model were sent to Roberts and Komesaroff, 
who used them in interpreting their recent observa- 
tions of Jupiter's decimeter radiation [Roberts and 
Komesaroff, 1964, 1965]. 

The purpose of this paper is three fold: (1) To present 
the numerical results which were used by Roberts 
and Komesaroff in interpreting their observations; 
(2) to present results of additional numerical calcula- 
tions based on the same model; (3) to use these new 
calculations to extend the interpretation begun by 
Roberts and Komesaroff. 

In addition to the results reported here, theoretical 
calculations of synchrotron radiation from Jupiter 
have been made by Chang and Davis [1962]; Korchak 
[1963], and by Ortwein, Chang, and Davis [1965]. 

2. The Model 

The model on which the calculations are based is as 
follows: (For additional details see [Thorne, 1963].) 



The magnetic field of Jupiter is idealized as a dipole 
field centered at the center of the planet. The radiat- 
ing electrons are confined to a narrow dipolar shell of 
equatorial radius equal to 3 Jovian radii and of thick- 
ness A/?. Shadowing by the planet is precisely con- 
fined to the region of the dipole shell hidden from the 
Earth by the planet's disk. At the equator the dis- 
tribution of electron energies, £*, and of pitch angles 
(angle between velocity vector and magnetic field), 
a ( > n is taken as 

A E-y 



N y (E, 0Le) = j£ —^r y 2 Of sin* <*e, 



if 



\E\ < E < E 2 , and 



iV v (£, Oe) = 0, if 



|a C rit< OL< (77 — CXcrit); 

(E<E U or E>E 2 , 
or a e < a CTiu or a e > (tt — a crit ); 



(1) 



(2) 



where Aq, y, E u E 2 



and a crit are constants. The 



constant a crit is chosen so as to rule out electrons with 
mirror points inside the planet. It is that pitch angle 
which corresponds to electrons with mirror points at 
Jupiter's surface 

a crit = sin- 1 (3- 7 / 4 ) = 8.4° 

[cf Thorne, 1963 (12) and (lb)]. The constant A fixes 
the total number of radiating electrons, and hence the 
intensity of the radiation. The constants E u £2, and 
y determine the electron energy distribution, and hence 
the spectrum of the radiation. The amplitude factors 
a q fix the electron pitch angle distribution, and hence 
the beaming and polarization of the radiation. 

This model differs from that of Thorne [1963] only 
in this, that shadowing by the planet's disk and the 
absence of electrons with mirror points inside the 
planet are taken into account. Comparison of nu- 
merical computations based on this model with com- 
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putations based on the original model reveals that 
shadowing and high mirror-point emission have only 
minor effects on the integrated intensity and polari- 
zation of the radiation. For instance, the total effect 
on the intensity is < 7 percent for ^ m in < ^ < ^max 
(cf©). 1 

3. Numerical Results 

Roberts and Komesaroff [1965] have confirmed the 
rinding of Gower [1963] that the spectrum of the de- 
cimeter radiation from Jupiter is flat 



h = dl/dv « p°. 



(4) 



If the decimeter region lies in the frequency range 
^min<^Omax, where 

^min=1.61 X 10 8 X 0B Ogauss ) X (£ lMeV ) 2 sec" 1 , 
*W = 1 -61 X 10 5 X (B 0e ) X (£ 2MeV ) 2 sec" 1 , (5) 



trons of Jupiter's Van Allen belt must be y=l. In 
table 1 we present data which enable one to calculate 
the intensity and polarization of the radiation as a 
function of observer's polar angle for y=l and 
^min <v< v max . I n particular, we give the intensity, 
Ivq, and polarization, P Q , for the special distributions 
of pitch-angle, 



CLq' bqq' 



1 if q' 4= q 
[1 if q' = q. 



(7) 



For the more general distribution (1) the intensity and 
polarization are (superposition principle) 



*v 2^ CLqlvq, 
Q 

P = (!//„) (2 W*). 



(8) 



and Bo is the equatorial field strength at 3 Jovian radii, 
then our model predicts [cf Thorne, 1963] 

/„ OC j,-(7-l)/2 ? ( 6 ) 

and consequently, the energy exponent for the elec- 



1 In figures 6 and 7 of [Thorne, 1963] the author presented estimates of the effects of 
shadowing and of high-latitude emission. It was not stated there, but it should have been, 
that those estimates are merely upper bounds and that for r e0 IR < 4 they are considerably 
larger than the true effects of shadowing and high latitude emission. Hence, figures 6 and 7 
are not useful in studying Jupiter's decimeter radiation, for which r e o/R ~ 3. 



The data of table 1 are not the data used by Roberts 
and Komesaroff [1965]. Since data for y= 1 were not 
available at the time of their analysis, they used the 
results of computations for y = 5/3, which are pre- 
sented in table 2. 

We do not know a priori that the decimeter region 
lies in the frequency range ^ min < v < max . Conse- 
quently, it is of interest to examine the beaming and 
polarization for the regions v < v min and v > v max . 
This is done in tables 3 and 4 for 7 = 5/3 (the analysis 
for other values of y is much more difficult) and for 
g = 3.3. 



Table 1. Intensity and polarization as functions of q and So for y = 1.0, v m in <v< v ma x a 







00 = 


90° 


6>o = 


87° 


0o = 83 o 


00 = 


80° 


00 = 


77° 




l vq 


P, 


h q 


P q 


/„„ 


P, 


h q 


P* 


h Q 


Pq 


1 


46.9 


-0.063 


47.1 


-0.062 


47.3 


-0.059 


47.4 


-0.056 


47.9 


-0.047 


2 


27.4 


+ .103 


27.4 


+ .103 


27.4 


+ .104 


27.4 


+ .103 


27.4 


+ .104 


3 


19.8 


.228 


19.8 


.228 


19.7 


.227 


19.6 


.225 


19.4 


.223 


4 


16.0 


.310 


15.9 


.310 


15.8 


.309 


15.7 


.307 


15.4 


.304 


5 


13.7 


.365 


13.7 


.365 


13.5 


.364 


13.3 


.361 


13.0 


.358 


6 


12.2 


.404 


12.1 


.403 


12.0 


.402 


11.7 


.399 


11.4 


.397 


8 


10.2 


.453 


10.1 


.452 


9.91 


.451 


9.60 


.448 


9.22 


.445 


10 


8.93 


.482 


8.87 


.482 


8.61 


.480 


8.25 


.478 


7.82 


.475 


13 


7.69 


.510 


7.62 


.509 


7.31 


.508 


6.90 


.506 


6.42 


.503 


16 


6.86 


.527 


6.77 


.527 


6.42 


.525 


5.97 


.523 


5.45 


.520 


20 


6.07 


.542 


5.97 


.541 


5.58 


.540 


5.09 


.538 


4.53 


.535 


25 


5.39 


.553 


5.28 


.553 


4.83 


.552 


4.30 


.549 


3.72 


.547 


30 


4.89 


.561 


4.77 


.561 


4.28 


.559 


3.72 


.557 


3.12 


.555 


50 


3.74 


+ .577 


3.59 


+ .577 


2.98 


+ .575 


2.36 


+ .573 


1.77 


+ .571 



a The quantities q and y are defined in (2) and (1); y m in and i> max are defined in (5); O is the angle between the mag- 
netic dipole axis and the observer's line of sight ("observer's polar angle") and l vq and P q are the intensity and polari- 
zation of the radiation for the special pitch-angle distribution (7). I vq is measured in units of 



k' 



■TT**^™- 



)(^-)( 



(MKS). 



where R is the radius of Jupiter and B is the equatorial field strength in the emission region (3 Jovian radii). This 
table was calculated from (26) and (27) of Thorne [1963] with obvious modifications for shadowing effects and for the 
absence of high mirror-point emission. The entries in the table are accurate to ± 1 percent. 
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Table 2. Intensity and polarization as functions of q and $ for by Roberts and Komesaroff [1965] are these: 

2.2XlO- 4 sin 3 a e cm- 3 e mA 
N v (E 9 a e )= u XL > ^~~MeV' Cm ' ^ 



y 


= 5/3 and v min < 


V < ^max 


a 






6»o = 90° 


0o = 


77° 


7 








Un 


Pfl 


/,,, 


P, 


2 


16.4 


+ 0.096 


16.5 


1 0.131 


3 


10.8 


.180 


10.6 


.174 


3.3 


9.86 


.220 


9.66 


.213 


4 


8.33 


.297 


8.06 


.288 


5 


6.99 


.374 


6.64 


.365 


6 


6.12 


.426 


5.73 


.417 


8 


5.06 


.490 


4.58 


.481 


10 


4.41 


.528 


3.86 


.519 


13 


3.78 


.562 


3.15 


.554 


16 


3.35 


.583 


2.66 


.574 


20 


2.96 


.600 


2.21 


.592 


25 


2.62 


.614 


1.81 


.606 


30 


2.38 


.623 


1.52 


.616 


50 


1.82 


+ .641 


0.856 


+ .634 



a Data used by Roberts and Komesaroff [1965]. Symbols 
footnote to table 1, except that l vq is measured in units of 



K J2!* &*(*$* (W „-„.(MKS). 

The entries in the table are accurate to ±1 percent. 



TABLE 3. Intensity and polarization as junctions of v\v m \ n and 0o 
for 7=5/3, ^r = 3.3, and v < v m ax ("Takeoff Region" ) a 





0,,= 


90° 


00 = 


77° 














/,„ 


/'</ 


/,„ 


1\ 


0.001 


0.47 


0.26 


+ 0.46 


+ 0.26 


.003 


.96 


.26 


.94 


.26 


.01 


2.06 


.268 


2.02 


.265 


. .03 


3.85 


.276 


3.7C) 


.273 


.1 


6.83 


.289 


6.67 


.284 


.3 


9.02 


.269 


8.81 


.262 


1.0 


9.71 


.232 


9.46 


.222 




9.86 


4 .220 


9.66 


+ .213 



a For explanation of symbols see footnot 


e to tab 


e2. 


Accuracy: ± 1 


per 


cent. Note that 


the entries in this table are related to t 


hose of 


tab 


e 4 by I(vlv min 


)tahl, 


.•* = [AW*w = 0) 


~/(WlW = 0.001 W^min)]table 4. 













TABLE 4. Intensity and polarization as functions of v\v mnx and O 
for y = 5/3, q = 3.3, and v > v m \n {"Takeoff Region")* 





00 = 


90° 


0„ = 


77° 














Iv* 


P* 


/,„ 


Pi 


0. 


9.86 


+ 0.220 


9.66 


+ 0.213 


1. 


9.39 


.218 


9.20 


.211 


3. 


8.90 


.216 


8.72 


.209 


10. 


7.80 


.208 


7.64 


.200 


30. 


6.01 


.184 


5.90 


.175 


100. 


3.03 


+ .065 


2.99 


+ .056 


300. 


0.835 


-.304 


0.852 


-.293 


1000. 


.151 


- .546 


.206 


-.196 



"For explanation of symbols see foot 



able 2. Accuracy: ±1 percent. 



4. Fitting of the Model to Jupiter 

The simplest distributions of electron energies and 
equatorial pitch angles which seem capable of repro- 
ducing the decimeter radiation of Jupiter as measured 



N v (E,Oe) = 



1.2X1 ()- 4 (sin 2 (x e + 2.0 sin 40 a e ) 



«„A« 



with ^ min < J> de cimeter< *W 

EWo £ 2 



E 

for 21 cm data; (9b) 
.e., with 
1() 4 , EW»<2. (10) 



Here AR is the thickness of the radiating shell as meas- 
ured in units of Jovian radii; B {) is the equatorial mag- 
netic field strength at 3 Jovian radii (the emission 
region) as measured in gauss; and the energies E, 
Eu and E z are measured in MeV. These electron 
distributions are capable of reproducing (1) the in- 
tensity, 6.7XK)- 26 Wmr 2 (c/s) -1 , of the radiation 
observed at Earth; (2) the spectrum of the radiation, 
l v a v°; (3) the beaming of the radiation, 

l v (13° magnetic latitude)//^ (magnetic equator) 

[ 0.89 at 21 cm 

- 1.0 at 74 cm; 

and (4) the degree of polarization of the radiation, 

P (magnetic equator) ~ 0.22 

(0.22 at 74 cm 
P (13° magnetic latitude) * \ (12) 

[0.18 at 21 cm. 

However, so long as the magnetic field is assumed to 
be that of a centered dipole, no electron distribution 
can account for the asymmetries in the radiation ob- 
served by Roberts and Komesaroff. 

Clearly, the electron distribution cannot be of the 
forms (9a) and (9b) simultaneously. Rather, as was 
first noted by Roberts and Komesaroff [1965J, the true 
distribution may be a mixture of these, in which the 
lower energy electrons have pitch angles similar to 
(9a) while the higher energy electrons have pitch angles 
similar to (9b). 2 

The electron distributions (9) are not completely 
unique; there are other distributions which will also 
reproduce the observed features of the radiation. 
However, the following limits can be put on the dis- 
tribution: 

Pitch angles. Any one-component model which 
fits the 74 cm data must have q ~ 3. For the 21 cm 
data, as was first noted by Roberts and Komesaroff 
[1965], one is forced to invoke at least two groups of 



-J. A. Roberts reported at this conference that the apparent dependence of pitch an<de 
on energy is in considerable doubt. Because his recent 48 cm data resemble closely the 
data at 21 cm. Roberts is inclined to disbelieve the low degree of equatorial beaming which 
seems to characterize his 74 cm data. 
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electrons — one group with pitch angles distributed 
nearly isotropically (q < 2) and the other with very 
flat pitches {q >20). 3 

Energy. The direction of polarization of Jupiter's 
decimeter radiation is observed to be orthogonal to the 
projection of the magnetic axis on the plane of the sky; 
but for v > 100 v max , our model demands (table 4; also 
Thorne [1963]) that the radiation be polarized parallel 
to the magnetic axis. Consequently, we can be cer- 
tain that decimeter ^ 100 v max , or, equivalently, that 
the high-energy cutoff in the electron energy distribu- 
tion satisfies 

[(£2) Me vPX(fio) g auss>200. (13) 

The more stringent bound, ^decimeter < ^max (E£Bo <£ 2 
X 10 4 ), seems quite likely, since the spectrum could be 
flat over a section of the region v > v max only if y < 1; 
and y < 1 seems unlikely on physical grounds. 

Although we are fairly certain that ^decimeter< ^max^ 
we cannot be sure that v min < ^decimeter ([^Mev] 2 #o gauss 
< 2). Table 3 reveals that the beaming and polariza- 
tion of the radiation in the region v < v min is not too 
different from that for the region ^ m i n < v < v max - 
However, the spectrum is quite different in these 
regions. If ^decimeter< ^mim then the energy exponent, 
y, must exceed 1 in order for the decimeter spectrum 
to be flat. In fact, y might be as large as 5 (the value 
for the Earth's Van Allen belt — [cf O'Brien et al., 
1962] — if ^ min were sufficiently large. Such a situation 
cannot be ruled out. 

s In a private communication Roberts and Komesaroff note that the particular 2-compo- 
nent pitch angle distribution given in their paper [Roberts and Komesaroff, 1965] to fit the 
21 cm data is not correct: but that the necessity for two components, one with q <2 and 
the other with q > 20. is unavoidable. 



The author is indebted to J. A. Roberts and M. M. 
Komesaroff for making the results of their observations 
available to him before publication. The numerical 
computations reported here were performed on the 
Princeton University IBM 7094 computer, which is 
supported in part by National Science Foundation 
prant NSF-GP579. The author was the recipient 
of fellowship support from the National Science Foun- 
dation, from the Danforth Foundation, and from the 
Woodrow Wilson Foundation during the period of 
this research. 
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Observations of Jupiter at 8.6 mm 

John E. Gibson 

E. O. Hulburt Center for Space Research, U.S. Naval Research Laboratory, Washington, D.C. 

Attempts to obtain measurements of the brightness temperature of Jupiter at 8.6 mm were 
made on several occasions near the opposition in May, 1959. With the 10-ft reflector the expected 
antenna temperature was low, and it was necessary to average repeated drift scans to obtain a 
measurable deflection. Atmospheric fluctuations nullified the results on some occasions, but 
analytical criteria found effective in more recent work have enabled the measurements for three 
days to be evaluated with some confidence. These results were 

1-2 May, 1959 308 ± 88 °K (p.e.) 

6-7 May 291±88°K 

8-9 June 260±90°K 

These values, obtained with north-south polarization, exceed the expected temperature by 
roughly a factor of two, and seem to indicate an anomalous effect in this period. The stated un- 
certainties were derived from the random fluctuations in the drift scans, with a relatively small 
systematic error not included. However, it should be emphasized that some chance exists that 
these results are in large error. For example, on another day, June 10, data considered to be com- 
parable in quality yielded a blurred unmeasurable result. Therefore, further observations at 

various periods would be useful. 

(Paper 69D12-590) 
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